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Abstract 



Both experimental and numerical studies of fluid motion indicate that initially 
localized regions of vorticity tend to evolve into isolated vortices and that these 
vortices then serve as organizing centers for the flow. In this paper we prove that 
in two dimensions localized regions of vorticity do evolve toward a vortex. More 
precisely we prove that any solution of the two-dimensional Navier-Stokes equation 
whose initial vorticity distribution is integrable converges to an explicit self-similar 
solution called "Oseen's vortex" . This implies that the Oseen vortices are dynami- 
cally stable for all values of the circulation Reynolds number, and our approach also 
shows that these vortices are the only solutions of the two-dimensional Navier-Stokes 
equation with a Dirac mass as initial vorticity. Finally, under slightly stronger as- 
sumptions on the vorticity distribution, we also give precise estimates on the rate 
of convergence toward the vortex. 

1 Introduction 

In this paper we consider the motion of an incompressible, viscous fluid in two-dimensional 
Euclidean space. The velocity of such a fluid is described by the Navier-Stokes equations 



where u = u(x, t) £ R 2 is the velocity field, p = p(x, t) £ R is the pressure field, and 
x £ R 2 , t > 0. For simplicity, the kinematic viscosity has been rescaled to 1. 

We prove two basic results about the solutions of (JTJ. First we show that for any 
initial velocity field whose vorticity is integrable, the solution of (JTJ with this initial 
velocity approaches an Oseen vortex, an explicit solution of (J^l exhibited below. As we 
also show, the Oseen vortex is in fact the unique solution of (JTJ) with a Dirac mass as initial 
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vorticity. We then examine in more detail the approach toward the vortex by studying 
the spectrum of the linearized equation around the vortex solution. If we assume that the 
initial vorticity distribution lies in a weighted L 2 space, we can derive estimates on the 
spectrum of this linearized operator which allow us to prove optimal bounds on the rate 
of convergence toward the vortex. 

We now describe our results in more detail. As we have argued in and it is 
often easier to understand the asymptotics of solutions of (JTJ) by studying the evolution of 
the vorticity, rather than the velocity. This is especially true in two- dimensions where the 
vorticity is a scalar. One can always then reconstruct the velocity field via the Biot-Savart 
law. Taking the curl of (JTJ we find that the vorticity u> = d\U 2 — d 2 U\ satisfies: 

^ + (u • V)w = Ac; , x G R 2 , t > . (2) 
The velocity field u is defined in terms of the vorticity via the Biot-Savart law 

u (*) = 7T I J—^ (V) x G R2 ' ( 3 ) 

2vry R 2 \x-y\ 2 

Here and in the sequel, if x = (xi, x 2 ) G R 2 , we denote x = (xi, x 2 ) T and x^ = (— x 2 , Xi) T . 

The vorticity equation is globally well-posed in the space L X (R 2 ). In particular, the 
results of Ben-Artzi j2], Brezis j7] and Kato [21] imply that: 

Theorem 1.1 For all initial data uo G L 1 (R 2 ) ; equation has a unique global solution 
lu G C°([0,oo),L 1 (R 2 )) n C ((0,oo),L oo (R 2 )) such that u(0) = u . Moreover, for all 
p G [1, +oo] ; there exists C p > such that 



K-,t)| P <^i, t>0. (4) 



C p \lUq\ 
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Here and in the remainder of the paper | • \ p denotes the norm on L P (R 2 ). If u G L q (R 2 ) 2 , 
we set \u\ q = | |u| \ q , where |u| = [u\ + u^) 1 ^ 2 . 

Among its other properties the semi-flow defined by (J2J) in L 1 (R 2 ) preserves mass, i.e. 



lu(x, t) dx = / u>q(x) dx , t>0. (5) 
r 2 Jn 2 

Furthermore, if the solution is sufficiently spatially localized so that the first moments of 
the vorticity distribution are finite then these are also preserved: 



XjOj(x, t) dx = / XjUi (x) dx , t > , j = 1, 2 . (6) 

R 2 </r 2 

It is important to realize that the solutions of (J2J) given by Theorem 11.11 correspond 
to infinite energy solutions of the Navier-Stokes equations ([!}. More precisely, if u(x,t) 
is a solution of (J2J) such that J uj(x,t) dx ^ 0, then the velocity field u(x,t) given by (jSJ) 
satisfies |u(-,t)| 2 = oo for all t. Explicit examples of such infinite energy solutions are the 
so-called Oseen vortices: 



u(x, t) = ^ e-IT/W , U (x, t ) = ~[l- e-W/m ) , (7) 
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where I +x 2 and a G R is a parameter which is often referred to as the "circulation 

Reynolds number". These solutions are "trivial" in the sense that u(x,t) ■ Vcu(x,t) = 0, 
so that (J2J) reduces to the linear heat equation. However, they play a prominent role in 
the long-time asymptotics of (J2J). Indeed, let 

G( {) = JL e -ia*, »o (£ ) = ^ii(i_ e -m«), 5eR =. ( 8 ) 

The following is the main result of this paper: 

Theorem 1.2 If Uq G L 1 (R 2 ), the solution u(x,t) of HJ) satisfies 



lim t 1 p 



a: 



u(;t) - jG(-j=) =0, for l<p<oo, (9) 



where a = f n2 uo(x)dx. If u(x,t) is the solution of (QJ) obtained from u(x,t) via the 
Biot-Savart law (OJ), then 



lim t' 2 i 



CY si 

u(-,t) - —= v G (— ) =0, for2<g<oo. (10) 



Vt 



In other words, the solutions of (J2J) in L X (R 2 ) behave asymptotically as the solutions 
of the linear heat equation d t uj = Au with the same initial data. For small solutions, this 
result has been obtained by Giga and Kambe in [T7], see also jEj. As was observed by 
Carpio [TU] (see also [IE]), there is a deep connection between the asymptotics Q and the 
uniqueness of the fundamental solution of the vorticity equation. More precisely, Carpio 
proved that © holds provided the Oseen vortex (J2J) is the unique solution of (j2J) with 
initial data ad, where 5 is Dirac's measure. As is shown in ^2], fU\ . this is true at 
least if | a | is sufficiently small. In this paper, we use a different method which allows us 
to obtain © without any restriction on a. As a by-product of our analysis, we prove the 
uniqueness of the solution of (0) with a (large) Dirac mass as initial condition: 

Proposition 1.3 Let u G C°((0, T), L\R 2 ) n L°°(R 2 )) be a solution of (0) which is 
bounded in L l (R 2 ), and assume that oo(-,t) (considered as a finite Radon measure on H 2 ) 
converges weakly to a5 as t —> 0+, for some a G R. Then 

ry / or \ 

w(x,t) = -G^J , xGR 2 , 0<t<T. 



Theorem ll.2l has a number of important consequences. Recalling that a can be thought 
of as the Reynolds number of the flow, Theorem 11.21 says in more physical terms that 
the Oseen vortices are globally stable for any value of this number. In contrast to many 
situations in hydrodynamics, such as the Poiseuille or the Taylor- Couette flows, increasing 
the Reynolds number does not produce any instability. From another point of view, our 
result is compatible with the "inverse cascade" of energy in two-dimensional turbulence 
theory. In contrast to the situation in three dimensions where energy injected into the 
system at large scales flows to smaller and smaller scales until it is dissipated by viscosity, 
in two dimensions both experimental and numerical results indicate that even for very 
turbulent, high Reynolds number flows, there is a tendency for smaller vortices to coalesce 
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and form larger and larger vortices. In this context, Theorem 11.21 savs that in the whole 
space R 2 this process continues until only a single vortex remains. 

Another consequence of this result is that the Oseen vortices are the only self-similar 
solutions of the Navier-Stokes equations in R 2 such that the vorticity field is integrable, 
see j2H] for a related result. For completeness, we mention that these equations have many 
other self-similar solutions with nonintegrable vorticities. Indeed, adapting the results of 
jH] to the two-dimensional case, it is easy to verify that the Cauchy problem for (JTJ is 
globally well-posed for small data in a Besov space which contains homogeneous functions 
of degree —1. For such initial data, the velocity u(x, t) and the vorticity ui(x,t) are 
automatically self-similar, due to scaling invariance. For instance, given any continuous 
map <f : S 1 — > R, there exists e > such that (JIJ has a self-similar solution u(x, t) with 
initial data u (x) = ex\x\~ 2 ip(x/\x\). If <p is not constant, the associated vorticity u(x,t) 
decays like l/|x| 2 as \x\ — > oo (at least in some directions), so that u>(-,t) £ L 1 (R 2 ). 

To prove Theorem 11.21 our strategy is to study a rescaled version of (J2J) which is sug- 
gested by the form of the vortex solution (J7j). Thus, we introduce the "scaling variables" 
or "similarity variables" : 

f=iL, r = log*. 

If u)(x,t) is a solution of (J2J) and if u(x,t) is the corresponding velocity field, we define 
new functions to(£, r), v(£, r) by 

u(x,t) = - w(Jj=,logf) , u(x,t) = -^v(-^,logt) . (11) 

Then r) satisfies the equation 

d T w + (v ■ V^)w = Cw , (12) 

where 

Cw = A%w+ -(£ • V()w + w . (13) 

The rescaled velocity v is reconstructed from the rescaled vorticity w by the Biot-Savart 
law: 

vK)= ill^¥'" ,( " )d "' « eR2 ' (i4) 

From (fTTj) and Theorem ll.l[ it is clear that the Cauchy problem for ()12j) is globally well- 
posed in L X (R 2 ). Remark that w(-, 0) = 1), hence imposing initial conditions to (fT2*|) 
at time r = corresponds to imposing initial conditions to (j2J) at time t = 1. This is of 
course harmless since (J2J) is autonomous. Observe also that the Oseen vortices {aG} ae R 
are by construction a family of equilibria of ()12l) . 

To prove Theorem 11.21 we first study the long-time asymptotics of solutions whose 
vorticity distribution is more strongly localized than is necessary just to be in L 1 . For 
any m > 0, we introduce the weighted Lebesgue space L 2 (m) defined by 

L\m) = {/eL 2 (R 2 )|||/|| m <oo}, where 

1/2 

= \b m f\ 2 ■ 



[ (i + l£lTl/(OI 2 d£ 

JR 2 
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Here and in the sequel, we denote &(£) = (1 + |^| 2 )^ 2 . If m > 1, then L 2 (m) L X (R 2 ). 
In this case, we denote by L^m) the closed subspace of L 2 {m) given by 

L 2 {m) = ju> G L 2 {m) J w (£)d£ = oj. (15) 

For v G (L 2 (m)) 2 , we set ||v|| m = || |v| || m , where |v| = (vf H-fl) 1 / 2 . 

As we observed in ^3] , a crucial advantage of using the vorticity formulation of the 
Navier-Stokes equations is that the spatial decay of the vorticity field is preserved under 
the evolution of (J2J. This remark is especially useful if one is interested in the long-time 
asymptotics since for parabolic equations there is a close relationship between the spatial 
and temporal decay of the solutions. The following result shows that the Cauchy problem 
for (fT2"|) is globally well-posed in the weighted space L 2 {m) if m > 1. 

Theorem 1.4 ( [14J , Theorem 3.2) Suppose that w$ G L 2 (m) for some m > 1. Then 
fflfy has a unique global solution w G C°([0, oo), L 2 {m)) with w(0) = w , and there exists 
C\ = Ci(||iUo||m) > such that 

\\w(r)\\ m < d , r>0. (16) 

Moreover, Ci(||iUo||m) ~ * as \\ w o\\m ~~ > 0. Finally, ifwo G L^m), then J R2 w(£, r) d£ = 
for all r > 0, and lim ||w(r)|| m = 0. 

Note that in contrast to this result the semi-flow defined by the Navier-Stokes equation 
does not preserve the spatial localization of the velocity field jTJ]. For instance, if the 
initial velocity u lies in L 2 (m) for some m > 2, then in general the solution u(-,£) of (|T|) 
with initial data u will not be in L 2 (m) for t > 0. For a detailed study of the localization 
properties of solutions of the Navier-Stokes and vorticity equations, we refer to the recent 
work of Brandolese [3] , jl] , [H] . 

Our proof of Theorem 11.21 begins with a proof that the Oseen vortices attract all 
solutions of (|12p with initial data in L 



m) 



Proposition 1.5 Let m > 1, wq G L 2 {m), and let w G C°([0, +oo), L 2 {m)) be the solu- 
tion of KWi) with initial data w . Then 

\\w(t) — aG\\ m — > as t —> +oo , 

where a = J R2 w (£) d^. 

Thus, any solution of the vorticity equation which is sufficiently localized to be in 
L 2 {m) for m > 1 will converge, as time tends to infinity, toward one of the Oseen vortices, 
regardless of how large the Reynolds number is. As is shown in Section El this global 
convergence result can then be extended to all solutions in L 1 (R 2 ). Returning to the 
original variables, we thus obtain Theorem 11.21 as a corollary. 

The proof of Proposition ll.5l is based on the existence of a pair of Lyapunov functionals 
for the rescaled vorticity equation (|12|). The first Lyapunov function, which is just the 
L 1 norm, is nonincreasing due to the maximum principle. It implies that the u;-limit set 
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of any solution must lie in the subset of solutions which are either everywhere positive or 
everywhere negative. 

On the subset of positive solutions, we use a second Lyapunov function which is 
motivated by a formal analogy between (|T2"j) and some kinetic models such as the Vlasov- 
Fokker-Planck equation. Given w : R 2 — > R + , we set 

^w = / Ri ^).og(||)d f . 

This quantity is often called the relative entropy (or relative information) of the vorticity 
distribution w with respect to the Gaussian G, see e.g. [30J. A direct calculation shows 
that H is non-increasing along the trajectories of (fT2"j) : 



—H(w) = — I 
dr 



R 



T-, , ( W 



2 



d£<0 



More precisely, this formula shows that H is strictly decreasing except along the family of 
Oseen vortices, which play the role of the Maxwellian distributions in kinetic theory. By 
LaSalle's invariance principle, the cu-limit set of any (nonnegative) solution of (|12j) with 
total vorticity a is thus reduced to a single point {aG}, which proves Proposition 11.51 

We next investigate the rate at which the solution w(t) of (|12j) approaches the vortex 
aG as r — > oo. This can be done by studying the linearization of ()12|) at the vortex. 
Under the assumptions of Proposition II .5[ there exists /i > such that — aG|| m = 

0(e~^ T ) as r — > oo. Convergence is thus exponential in the rescaled time r = \ogt, 
hence algebraic in the original time t. As in the case of the linear heat equation, the 
convergence rate satisfies fi < (m— 1)/2. This limitation originates from the essential 
spectrum of the linearized operator at the vortex, and is related to the spatial decay 
of the vorticity. In addition, discrete eigenvalues prevent convergence to be arbitrarily 
fast even for exponentially localized solutions. The following statement generalizes to 
arbitrary data in L 2 (m) the local results of [T3| . 

Proposition 1.6 Fix m > 2. For any wq G L 2 (m), the solution w G C°([0, +oo), L 2 (m)) 
°f lIBj) with initial data wo satisfies 

\\w(t) -aG\\ m = 0(e- r/2 ) , as r -> oo , (17) 

where a = J R2 wo(£) d£- Moreover, if m > 3 and (3\ — fa — where fa = J R2 £iWo(0 d£, 
then 

\\w(t) — aG\\ m = 0(e~ T ) , as t — >• oo . (18) 

As we shall see in Section|H if m > 2 and (fa, fa) ^ (0, 0), then ||w(r) — «G|| m decays 
exactly like e _T//2 as r — > oo, so that p7|) is sharp. Similarly, if m > 3, f)18|) is sharp in the 
sense that there is in general a correction to the Gaussian asymptotics decaying exactly 
like e~ T as r — > oo. If a 7^ 0, there is no loss of generality in assuming that fa — fa — 0, 
since this can always been achieved by an appropriate choice of the origin in the original 
variable x G R 2 (see Section HJ. The situation is different if a = 0, see [Hj. In this case, if 
{fa, fa) 7^ (0,0), the solution converges to zero at the rate e~ T//2 , and the next correction 
is Oije~ T ) due to secular terms in the asymptotics. 
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If we reexpress estimate (J17J) in terms of the original dependent and independent 
variables, we find that 



0(t 



u 



a 



r G( 



y/t \/i 



("7?) =°( t 



as t — - > oo, which represents a considerable sharpening of the decay rates in ©. Similarly, 



(118)) implies that the quantities above are 0(t ^ p^) and 0{t 



.(3. 
^2 



respectively. 



The rest of this paper is organized as follows. In Section |21 we study the compactness 
properties of the solutions of (fT2*j) in both L 2 (m) and L 1 (R 2 ). In Sectional we show 
that ()12|) has two Lyapunov functions, and we prove Proposition 11.51 Theorem 11.21 and 
Proposition 11.31 Finally, in Section HI we study the spectrum of the linearization of (|T2|) 
at the Oseen vortex and we obtain precise estimates of the rate at which solutions of (fT2j) 
approach the vortex. In particular, we prove Proposition 11.61 

Acknowledgements. The first author is indebted to J. Dolbeault and, especially, to 
C. Villani for suggesting the beautiful idea of using the Boltzmann entropy functional in 
the context of the two-dimensional Navier-Stokes equation. The research of C.E.W. is 
supported in part by the NSF under grant DMS-0103915. 



2 Smoothing and compactness properties 

In this section we discuss some general properties of solutions of the vorticity equation 
(|12|) which we will need to establish the convergence results of the Section El To control 
the nonlinear terms in (J2J) or (|12j) . we will need estimates on the velocity in terms of 
the vorticity. Let uo and u be related via the Biot-Savart law Q. If uo G L P (R 2 ) for 
some p G (1,2), it follows from the classical Hardy-Littlewood-Sobolev inequality that 
u G L q (~R 2 ) 2 where - = - — £, and there exists C > such that 

|u|, < C\u\ p . (19) 

Of course, a similar result holds if w and v are related via ([14)1 . Further estimates are 
collected in ([14). Lemma 2.1 and Appendix B). 



2.1 Compactness in L 2 (m) 

If Wo G L 2 (m) for some m > 1, we know from Theorem 11.41 that equation ()12|) has a 
unique global solution w G C°([0, oo), L 2 {m)) with initial data w . As is explained in [T4] . 
w is in fact a solution of the associated integral equation 



w (t) = S(t)w - S(t - s)v(s) • Vw(s)ds (20) 
Jo 

= S{t)w - I e -!( T -s) V -S{t- s)v{s)w{s)ds , 
Jo 

where S(r) = exp(r£) is the Co-semigroup generated by the operator £. Remark that, 
since S(t) is not an analytic semigroup, the solution w of ((2D|) is not (in general) a classical 
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solution of (|12|). In particular, r i— > u>(r) is not different iable in L 2 (m). For later use, we 
recall the following results ([E], Appendix A): 

1. If m > 1, there exists C > such that, for all / G L 2 (m), 

\\S(r)f\\ m < C\\f\\ m , \\VS(r)f\\ m <-^\\f\\ m , r>0, (21) 
where a(r) = 1 — e~ T . 

2. If < n < 1/2 and m > 1 + 2/i, there exists C > such that, for all / G Lg(m), 

\\S(r)f\\ m <Ce-^\\f\\ m , ||V5(r)/|U<^73^||/|| m , r>0. (22) 

3. Let 1 < q < p < oo and T > 0. For all a G N 2 there exists C > such that 

r<9 Q S(r)/| p < JL^^/j,, 0<r<T, (23) 

where 6(f) = (1 + |£| 2 ) 1/2 . 

From Theorem 11.41 we also know that there exists K\ > such that ||w(r)|| m < K\ for all 
r > 0. The aim of this section is to prove that the trajectory {u>(t)} t > is in fact relatively 
compact in L 2 (m). This is because equation ()12j) is both regularizing and "asymptotically 
confining", in the sense that solutions of (|T2"|) in L 2 {m) are asymptotically (as r — ► +oo) 
well localized in space. The localization effect originates in the dilation term |£ • V in the 
linear operator C, and hence does not occur in the original vorticity equation (j2J). 

We first prove that, for positive times, the solution w{r) belongs to the weighted 
Sobolev space 

H l {m) = {w G L 2 {m) \ diW G L 2 {m) for i = 1, 2} , (24) 
which we equip with the norm ||i/j||#i( m ) = (HHIm + ll^^llm) 1 ^ 2 - 

Lemma 2.1 Let wq G L 2 (m) with m > 1, and let w G C°([0, oo), L 2 (m)) be the solution 
of KWi) with initial data w . Then there exists K 2 > such that 

\\yw(r)\\ m <^- , forallrX), (25) 

where a(r) = 1 — e~ T . 

Proof: Consider the Banach space X = C°([0, T], L 2 (m)) n C°{{0,T], H^m)) equipped 
with the norm 

\\w\\ x = sup ||w(r)|| m + sup a(r) 1/2 ||Vw(r)|| m . 

re[0,T] re(0,T] 

We shall prove that there exist T > and K > such that, for all initial data Wq G L 2 {m) 
with 1 1 Wo || m < K\, equation (j2Tj|) has a unique solution w G X, which satisfies \\w\\x < -f^- 
We then apply this result to w = w(nT/2) with n G N. By uniqueness, we have 
?D(t) = w(t + nT/2) for r G [0,T], hence 

sup a{r) 1/2 \\Vw{T + nT/2)\\ m <K, for all n G N , 

0<t<T 
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which implies ()25p . 

To prove that (|2(J|) has a unique solution in X, we proceed as in ((Ej, Lemma 3.1). 
Given w%, w 2 & X, we define 

R(t)= ( 5(r-s)vi(s) • Vw 2 (s)ds , 0<r<T, 

where Vi(s) is the velocity field obtained from wi(s) via the Biot-Savart law (jl4j) . If 
q £ (1,2), we know from that 



It remains to estimate 



\\R(r)\\ m <c^ —L-^iKiuilt^iu 

VR(t)= [ V<S(r - s)vi(s) • Vw 2 (s)ds . 
Jo 



Applying (j23j) with p = 2 and g as above, we obtain 



|V<S(r - s) Vl (s) ■ Vw 2 (s)|U < , ^ u/J ^ViOOVwzOOIg 

a r — s 1/9 



Next, using Holder's inequality, estimate (|T^|) and the fact that L 2 (m) •—>■ L 9 (R 2 ) when 
2/(m+l) < g < 2, we find 

|& m Vl (s) • Vw 2 (s)| g < |v 1 (s)|^ L |6 m Vw; 2 (s)| 2 < CK(s)|,||Vw 2 (s)|| m (26) 

< C||wi(s)|| m ||Vw 2 (s)|| m . 

Thus 

a(T)^||Vi2(T)|| m < f f^Xj I^^llmllV^IUda 

y a(r-s) 1 /? 

* c (r ^-iK W i /- d ') K|uiwu - <2T) 

Summarizing, we have shown that ||-R||x < ^COI|wi||x||w 2 ||x, where C(T) — > as T — > 0. 
Using this bilinear estimate together with the bound (|21|) on the linear semigroup, it is 
straightforward to verify by a fixed point argument that, if T > is sufficiently small, 
then ()2()j) has a unique solution in X with the desired properties. □ 

Next, we decompose w(£,r) = qG(() + R(£,t), where G(£) is given by (JHJ) and a = 
Jr. 2 r ) d£ * s time- independent due to (J3J). Then R(-,t) G -^o( m ) f° r all r > 0, where 
Lg(m) is defined in ([T5]). Moreover, R satisfies the evolution equation 

d T R = CR- aAR - N(R) , (28) 

where 

A_R = v G • VR + v R ■ VG , N(R)=v R -VR. (29) 
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Here v G is given by (JKJ) and v R denotes the velocity field associated to the vorticity R by 
the Biot-Savart law (|14|) . The corresponding integral equation is 

R( T ) =S(T)R -a / S(t - s)AR(s) ds - / S(r - s)N(R(s))ds . (30) 
Jo Jo 

Equations (fl"2*j) and ()28|) are clearly equivalent. In particular, Theorem 11.41 implies that, 
given any a G R, the Cauchy problem for (|2"%|) is globally well-posed in L^{m) if m > 1. 

We now prove that positive trajectories of ()12j) in L 2 (m) are relatively compact. 

Lemma 2.2 Let w G L 2 (m) with m > 1, and let w G C° ( [0, oo), L 2 (m)) be the solution 
of Mty) with initial data Wo. Then the trajectory {w(r)} T > is relatively compact in L 2 (m). 

Proof: Since w G C°([0, oo), L 2 (m)), it is sufficient to show that {w(r)} r >i is relatively 
compact in L 2 {m). We decompose w{^,t) = aG(£) + R(£,t) as above, and consider 
the integral equation (jBllj) satisfied by R. Let -Ri(t) = R(t) — S(t)Ro, so that w(r) = 
aG + S(t)R q + Ri(t). Since R G Lp(m), it follows from that S(t)Rq converges 
to zero in H x (m) as r — > oo. In particular, since {w(r)} T > 1 is bounded in H 1 {m) by 
Lemma l2~Tl we see that {-Ri(r)} r >! is also bounded in H l {m). Now, we shall prove that 
{i?i(r)} T > is bounded in L 2 (m + 1). By Rellich's criterion (see [27], Theorem XIII. 65), 
this implies that {Ri(r)} T >i is relatively compact in L 2 {m). Since S(t)R q converges to 
zero as r — > oo, it follows that {w(t)} t >i is also relatively compact in L 2 (m), which is 
the desired result. 

To prove the claim, we decompose Ri(t) = R2ij) + R^ij) + R^ij) where 

R 2 ( T ) = -a [ e~^ T ~ s) V -S(T-s)v G R(s)ds , 
Jo 

R 3 ( T ) = - a I e~^ T - s) V -S{t - s)v R {s)Gds , 
Jo 

R 4 (t) = - [ e-^ r - s V ■ S(r - s)v R (s)R(s) ds . 
Jo 

To bound R2(t), we observe that 6v G G L°°(R 2 ). Using (|2"T]). we thus find 

\\MT)\\m+i = \a\ b m+1 [ T e~^- s ^-S(T-s)v G R(s) ds 

Jo 2 

< C\a\ f e-^ T - s) -\b m+1 v G R(s)\ 2 ds (31) 

Jo a(r — s) 2 

< C\a\ f e-^ T ~ 8 ) - 1 ^ |fo%||fl(a)|| m da<glff , 

Jo a(r — s) 2 

for all t > 0, where i£i = sup T>0 ||w(r)|| m . To bound i2a(r), we proceed in the same way. 
Using Holder's inequality, estimate (JT§j) and the embedding L 2 (m) °-> L 4 / 3 , we obtain 

|6 m+1 v R G| 2 < |v R | 4 |6 m+1 G| 4 < C\R\ A/3 < C\\R\\ m , (32) 
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so that ||i2 3 (r)|| m+ i < CK\ for all r > 0. Finally, to bound R±{j) we fix q G (1,2) such 
that q > 2/m. Combining (|21j) and (|23j) . we obtain 

\\Ra(t)\Ui < C £ e-ilr-O l/q \ b m+1 v R (s)R(s) \ q ds . (33) 

By Holder's inequality, \b\ R b m R\ q < |6v R |^2_ ||-R|| m . Applying Proposition B.l(2) in [TJ], 
we also obtain \bv R \2q_ < C||i2||2 < C||i2|L. Inserting these bounds into (JSHJ), we find 

2-9 9 

that ||i?4(r)|| m+ i < CKf for all r > 0. Summarizing, we have shown that ||-Ri(r)|| m+ i < 
CK\ for all r > 0, which concludes the proof. □ 

Finally, we show that negative or complete trajectories of (fT2"|) in L 2 (m) that are 
bounded in L 2 (m) for some m > 1 are also relatively compact. 

Lemma 2.3 Assume that m > 1 and that w G C°(R, L 2 (m)) is a solution of liij)) which 
is bounded in L 2 (m). Then {w(t)} t£ r, zs relatively compact in L 2 (m). 

Proof: By assumption, there exists K\ > such that ||iu(r)|| TO < i^i for all r G R. 
As in the proof of Lemma l2~2l we decompose w(C,,t) = aG((,) + i?(^,r), where a = 
Jr. 2 w (£> r ) d£- The remainder i?(r) satisfies the integral equation 

R( T ) = S(t - t )R(t ) -a S(r - s)AR(s) ds - S(t - s)N(R(s)) ds , (34) 

J TO J TO 

for all t < r. Since i?(r ) G -^o( m ) anc ^ ll-^( r o)l|m < CK 1 for all r G R, it follows from 
(f22|) that ||«S(r — r )i?(r ) || m — > as r — > — oo. Moreover, proceeding as in the proof of 
Lemma 12.21 and using the analogues of estimates (}3"Tj) , (|3*3*jl , it is easy to see that both 
integrals in the right-hand side of (|34|) have a limit in L 2 (m) (and even in L 2 (m + 1)) as 
r — > — oo. Thus, we have the representation 

/T t*T 
S{t- s)AR(s)ds- S(t - s)N(R(s))ds , tgR, (35) 
-oo J — oo 

which implies that ||i?(r)|| m+1 < CKf for all r G R. This shows that {w(r)} reR is 
bounded in L 2 (m + 1). On the other hand, it follows from Lemma f2. II that {w(r)} rG R is 
bounded in H l (m), hence {w(t)} t£ h is relatively compact in L 2 (m) by Rellich's criterion. 

□ 

Remark 2.4 By a bootstrap argument, it is clear from the proof of Lemma \2.tA that 
{w(r)} re R is bounded in H k (m') for all k G N and all m' G N. In other words, the 
trajectory {w(t)} t( zr is bounded in the Schwartz space S(R 2 ). 

2.2 Compactness in L X (R 2 ) 

We now study the compactness properties of the solutions of (|12|) in L X (R 2 ). We first 
recall two important estimates for the solutions of the original vorticity equation (J2J), and 
use the change of variables (|TT|) to obtain the corresponding bounds on the solutions of 

G2). 



n 



The first is a smoothing estimate, see for instance [21]. Under the assumptions of 
Theorem II .![ the vorticity u(x,t) satisfies for all p G (1, oo] 

C 

\Vu(;t)\ p <-aZT , t>0, (36) 
t 2 p 

where C' p depends only on |a>o|i. The second is a nice pointwise estimate due to Carlen 
and Loss, see Theorem 3 in [9.. For any (3 G (0, 1), there exists Cp > (depending on 
l^o | i) sucn that 

r i \x P 

\co(x,t)\<C p J -exp(-p lX ~ t yl )\u (y)\dy, x G R 2 , t>0. (37) 

Assume now that w Q G L^R 2 ), and let w G C°([0, oo), L : (R 2 )) be the solution of (fT2|) 
with initial data Wq. Then the function u(x,t) defined by (|TT|) is a solution of (J2J) on the 
time interval [1, oo) with initial data lu(x, 1) = Wo(x). Applying (JSfiJl and returning to the 
rescaled variables, we obtain 

|Vw(-,r)| p < °l ; , r>0, (38) 
a(r) 2 p 

where a(r) = 1 — e~ r . Similarly, we deduce from (|37|) that 

\w({,t)\<C p [ -Lexp(-[3 ^~ ye T/2l2 )\w (y)\dy , (6R 2 , r>0. (39) 
J R 2 a(r) V 4a(r) / 

With these estimates at hand, we are now ready to prove the main result of this section. 

Lemma 2.5 Let w G -^(R 2 ), and let w G C°([0, oo), L X (R 2 )) be the solution of KTty 
with initial data wo- Then {w(r)} T > is relatively compact in L : (R 2 ). 

Proof: Again, it is sufficient to show that {w(t)} t >i is relatively compact. Fix R > 0. 
From one has \w(£, r)\ < W\(£, r) + w 2 (C, r) where 

u>i(£,t) = Cp I — ^expf-ff ^ ye - ) \w (y)\ dy , 



|2/|<-R 



a(r) V 4a(r) 



w 2 (£,t) = Cp — -exp(-P )\w (y)\dy. 



4a(r) 



If |f | > 2R, then |f - ye~ T/2 \ > |£| - \y\ > |f |/2 whenever \y\ < R and r > 0. It follows 
that 2 

^ 1 (e,T)<^| Wo |i^ex P (-/3 T ^L-) , \£\>2R, 



a(r) V 16a(r). 
/ wi(f,r)df < Cp\w \i [ exp(-/?i£) dz =: e x (i2) . 

J\f\>2R J\z\>2R V 1,3 y 



hence 
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Moreover, using Fubini's theorem, we find 



Jr 2 P J\y\>R 



Thus we have shown 



sup/ \w(£,r)\d£ < e^R) + e 2 (R) -»• as R -> oo . (40) 

t>0 J\£\>2R 

On the other hand, using (f3"HJ) with p = oo, we see that 

K : = sup sup |Vu>(£, t)| < oo . 

Fix e > 0. According to pUj). there exists i? > 1 such that 

sup / H£,r) |d£ < e/3 . 
t>1 J\£\>R-1 

Let 5 = min(l, e(3KnR 2 )~ 1 ). If y e R 2 satisfies |y| < 5, then for all r > 1 



/ me-y,r)-w(e,r)id^ < 2 / me,T)id£< 

</|£|>.R -'lf|>-R-l 



2£ 

y 



and 



\w(Z-y,T)-w(Z,T)\d£ < 7rR 2 \y\ sup \Vw(£,t)\ < ttR a 5K < - , 
|?|<-R |5|<H+i d 

so that J R2 |w(£ — y,r) — w(£, t)| d£ < e. Thus we have shown 

sup sup / |u>(£ — y, t) — w(£, r)| d£ — > as 5 — > . (41) 

t>1 |y|<5jR2 

By the Riesz criterion (|27j. Theorem XIII. 66), it follows from (|4()jl . (f4*T]) that {w(t)} t >i 
is relatively compact in L X (R 2 ), which is the desired result. □ 

To conclude this subsection, we remark that Lemma 12.31 has no analogue in L X (R 2 ), 
namely negative trajectories of (fT2|) that are bounded in L X (R 2 ) need not be relatively 
compact in that space. To see this, let uo G L X (R 2 ) and let u(x,t) be the solution of (0) 
given by Theorem ll.il If r) = e r u;(£e r / 2 , e T ), then {w(-, r)} r < is a bounded negative 
trajectory in L 1 (R 2 ). However, w(-,r) is "evanescent" as r — > — oo, hence {w(-,t)} t < is 
not relatively compact in L X (R 2 ) unless co> = 0. 



2.3 Preservation of positivity 

A more qualitative property that will be essential for our analysis is the fact that solutions 
of (J2J), and hence of (|12|). satisfy a maximum principle. We state this property in the 
original variables, and for a generalized version of (J2J) where the velocity and vorticity 
field are not necessarily connected to each other. This generalization will be useful in the 
next section. 
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Proposition 2.6 Assume that u G C 6 °(R 2 x [0, oo), R 2 ) and that u G C 2 (R 2 x [0, oo), R) 
is a solution of 

d t u(x,t) + u(x,t) ■ Vu(x,t) = Au(x,t) , xeR 2 ,t>0. (42) 

lfu(x,0) > /or x G R 2 , then either u(x,t) = or u(x,t) > /or a// x G R 2 and a// 
t > 0. 

Proof: This classical result is obtained for instance by combining Theorems 3.5 and 3.10 
in the book of Protter and Weinberger • D 

As a corollary, we obtain that (0) preserves positivity. The same property holds for 
(fT2"j) and is proved using the change of variables (JTTJ). 

Corollary 2.7 Assume that ujq G L 1 (R 2 ) satisfies uq(x) > for almost all x G R 2 ; 
and that ojq{x) does not vanish almost everywhere. Then the solution of |1J) given by 
Theorem \l.l\ satisfies u(x, t) > for all x G R 2 and all t > 0. 

Proof: If u G <S(R 2 ) and to (x) > for all x G R 2 , the solution u>(x,t) of © and 
the corresponding velocity field u(x, t) satisfy the assumptions of Proposition 12.61 hence 
u)(x,t) > for all x G R 2 and all t > 0. Since solutions of J2J) depend continuously on 
the initial data in L 1 (R 2 ), the same result holds for any solution of (J2J with initial data 
ujq G L x (R 2 ) such that u (x) > almost everywhere. Moreover, given any to > 0, this 
solution satisfies uj G C 2 (R 2 x [t ,oo),R) and the corresponding velocity field satisfies 
u G C°(R 2 x [t , oo), R 2 ). Applying Proposition 12.61 again, we deduce that u(x, t) > for 
all igR 2 and all t > t , unless u(x, t ) = 0. Since this is true for any t > 0, we conclude 
that u(x,t) > for all x G R 2 and all t > 0, unless ojq{x) vanishes almost everywhere. □ 

3 Global convergence results 

This section is devoted to the proofs of Proposition II .51 Theorem 11.21 and Proposition II .31 
The argument relies on the compactness properties of the previous section, and uses the 
crucial fact that system (|T2*j) has two Lyapunov functions which we introduce now. 

3.1 A pair of Lyapunov functions 

Let $ : L X (R 2 ) — > R + be the continuous function defined by 

$H = f H0|d£, (43) 

and let 




In words, a function w G L : (R 2 ) belongs to S if and only if w(£) has (almost everywhere) 
a constant sign. Remark that, by Corollary 12.71 the set S is positively invariant under 
the evolution of (JT2J). 

We first show that $ is a Lyapunov function for the semiflow of (|12|) . More precisely, 
$ is strictly decreasing along the trajectories of (|T2"|) except on the invariant set £ where 
$ is constant. 
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Lemma 3.1 Let w G L^R 2 ), and let w G C°([0, oo), L X (R 2 )) be the solution of Iffify 
with initial data w . Then $(w(r)) < $(u>o) for all r > 0. Moreover $(w(r)) = §(w ) 
for all r > i/ and on/?/ i/wo G E. 

Proof: If w G E, then to(r) G E for all r > 0. Using (jSj), we thus find 



<&(w(t)) 



u>(£,r)d£ 



R 2 



R 2 



$Oo) 



for all r > 



Assume now that w £ E. Then u> = n>o~ — n> > where 

<(0 = maxKK), 0) > , w (O = - wm(w (£), 0) > . 

By assumption, both Wq and Wq are nonzero on a set of positive Lebesgue measure. Let 
wi and u>2 be solutions of 



d T wi + v ■ Vwi = Cw\ , t > 

9 T W 2 + V • V^2 = C-W2 , T > 



(44) 



with initial data u>i(0) = Wq~, w 2 (0) = W7 , where v(£, r) is the velocity field associated 
with the solution r) of (|12j) . Following the proof of Theorem ll .11 (see |21|). one verifies 
that iui,tt;2 G C°([0,cx)),L 1 (R 2 )) and that 

w 1 (z,T)dt= [ <(e)de, / W2 (e,r)de= / ^o(e)de, (45) 

r 2 Jn 2 Jr. 2 Jr. 2 

for all r > 0. Moreover, using Proposition 12.61 and a density argument as in the proof 
of Corollary 12.71 it is straightforward to verify that Wi(£, r) > and u> 2 (£,r) > for all 
(6R 2 and all r > 0. Now, by construction we have u>(£, r) = u>i(£, r) — u>2(£, r ), hence 

H£, r )l = l^i(^' r ) - w 2 (^,r)| < w 1 (£,t)+w 2 (£,t) , 
for all £ G R 2 and all r > 0. Integrating over R 2 and using (J43|) . we obtain 



/ K£,r)|d£ < / K(e,r)+ W2 (e,r))d£ 

JR 2 JR 2 



R 2 



K + (0 + ^ "(0)d£ 



R 2 



MOIdC, r>0. 



□ 



This shows that $(w(t)) < $(w ) for all r > 0. 

Next we fix m > 3 and we consider solutions of ()12j) in the invariant cone L 2 {m) n E + , 
where 

E + = {w G L 1 (R 2 ) | u>(£) > almost everywhere} . 
We define H : L 2 (m) n E+ -> R by 



. 



^(Oiog^K. 



(46) 
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Since 



wlog^j = w\og(ATTw) + w < Cw 2 + w , 

it is clear that H is well-defined and bounded from below by — 1/e. Moreover, using for 
instance the inequality 

\wi\ogwi - w 2 logw 2 \ < c(\w\ - w 2 \ 1/2 + \wi - w 2 \(w\ /2 + wl /2 )^j , 

one verifies that H is continuous on L 2 (m) H S + (equipped with the topology of L 2 (m)). 
We now show that is indeed a Lyapunov function for the semiflow defined by (|12|) on 
L 2 (m) n £+. 

Lemma 3.2 Assume thatwo G L 2 (m) flS + wrat/i m > 3, and let w G C°([0, oo), L 2 {m)) be 
the solution of Alty) with initial data wq. Then H(w(r)) < H(wq) for all r > 0. Moreover 
H(w(t)) = H(wq) for all r > if and only if wq = aG for some a > 0. 

Proof: If wq = aG for some a > 0, then w(r) = aG for all r > 0, hence obviously 
H(w(t)) = H(wq) for all r > 0. Assume now that wq is not a multiple of G (in particular, 
w 7^ 0). Then the solution ty(£, r) of (fT2J) is smooth and strictly positive for all r > 0. 
We claim that r t— > H(w(r)) is differentiable for r > 0, and that 

A#( w ( r )) = -J( w ( r )) , r > , (47) 
dr 



where 



IH=/ w(0 Vlog(^) \>0 . (48) 



Remark that I(w) vanishes if and only if w is proportional to G. Thus, under the as- 
sumptions above, it is clear that I(w(r)) > at least for r > sufficiently small, hence 
H(w(t)) < H(w ) for all r > 0. 

Thus, all that remains is to prove (|47p. Assume first that wq belongs to the Schwartz 
space S(R 2 ). Then w G C\[0, oo),S(R 2 )) is a classical solution of (JUJ) in S(R 2 ). More- 
over, w(£, t) satisfies a Gaussian lower boundfor any r > 0, see j21] or ([IE]) Theorem 3.1). 
More precisely, there exist positive constants 7 and C 1 (depending only on \wq\i) such 
that, for all £ G R 2 and all r > 0, 

w(£,r) > — - / expl-7 — )w {y)dy, 

a i T ) Jn? v 2a(r) / 



> -^T ex p(-7^)/ exp(- 7 J^W2/)dy, (49) 
a(r) V a(r)/J R2 V a[r)J 

see also (|3*9^l . Using these properties, it is straightforward to verify that r \— > H(w(t)) is 
differentiable for r > 0, and that 

^Z#Mr)) = J (l + log |) <9 r w d£ = jT ^ (l + log |) (£ W - v • Vw) d£ . 
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Next, using the identity Cw = div(GV(^)) and integrating by parts, we obtain 

/ R2 (i + io g |) ( ^de = -/ R2 v(io g |)-§v(|>d, 



/ R H v ( log I 



de = -i(w) 



On the other hand, using v • Vw = div(vw) and integrating by parts, we find 

/ (l + log^)(vVw)df = / (l + log(47™))(v-Vw)d£+ / ii!!(v.Vw)d£ 

JR 2 V ^ JR 2 JR 2 4 

= - / v- Vwd£-^ / {$-v)wd£. 

JR 2 1 </R 2 

We claim that both integrals in the last expression vanish. This is obvious for the first 
one, since v • Vw = div(viw). As for the second one, using (JHJ) and Fubini's theorem we 
obtain 



Jr 2 27r JR 2 xR 2 IC-^r 



Summarizing, we have shown that 



if Mn)) - F(w(r )) = - / I(w(r)) dr , (50) 



for all T\ > t > 0. 

We now return to the general case where w G L 2 {m) PI £+. Given T\ > r > 0, the 
solution w of (|12|) satisfies w G C°([r , Ti], H k {m)) for any G N, where H k {m) is the 
weighted Sobolev space defined in analogy with (|24|) . Moreover, the map 

w G L 2 (m) ^ we C°([r , n], ff fc (m)) 

is continuous. On the other hand, using for instance [22]; h is not difficult to verify 
that the quantity I(w) is finite for any positive w G H k {m) if k > 2, and that f(w) 
depends continuously on u> in that topology. Thus we see that both sides of (|50|) depend 
continuously on the initial data wo in the topology of L 2 {m). Since (|50|l holds for all wq 
in the dense subset 5(R 2 ) fl S + , it follows that (|5L)p is valid for all w G L 2 {m) R S + and 
all ri > r > 0. This concludes the proof. □ 

3.2 Convergence in L 2 (m) 

Using the compactness properties of Section |2] and the two Lyapunov functions of the 
previous subsection, we are now able to prove Proposition 11.51 

Lemma 3.3 Assume that m > 1 and that {u>(r)} rg R is a complete trajectory of Alty) 
which is bounded in L 2 {m). Then w{r) = aG for all r G R, where a = J R2 w(£, 0) d£. 
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Proof: We know from Lemma l2~3~l that {w(r)} Tg R is relatively compact in L 2 (m). In view 
of Remark 12.41 we can assume without loss of generality that m > 3. Let Q be the tu-limit 
set of the trajectory {if(r)} reR . Since by Lemma I3~T1 $ is a Lyapunov function which is 
strictly decreasing except on £, it follows from LaSalle's invariance principle that flc E. 
In particular, since the total mass is conserved, any w G Q satisfies $(iZ>) = | J R2 w(£) d£| = 
\a\. The same is true for any function w in the a-limit set A. As r i— > $(io(t)) is non- 
increasing, it follows that $(w(r)) = |a| for all r G R. By Lemma \3 . II again, we conclude 
that u>(t) G £ for all r G R. Thus, upon replacing w(£ 1 ,£ 2 ,' r ) by — w(£ 2 , £i> t), we can 
assume that {w(r)} reR C L 2 {m) R £+. 

We now use the second Lyapunov function if. By Lemma I3~2l and LaSalle's principle, 
A and f2 are contained in the line of equilibria {a'G} a '>o. Since the total mass is conserved, 
we necessarily have A = Q = {aG}. As H is non-increasing, it follows that H(w{r)) = 
H(aG) = alog(a) for all r G R. By Lemma I3.2I again, we conclude that w(t) = aG for 
all r G R. □ 

Proof of Proposition lTT^l Let w G L 2 {m) with m > 1, and let w G C°([0, +oo), L 2 {m)) 
be the solution of ()12|) with initial data Wq. From Lemma 12.21 we know that {u>(t)} t > 
is relatively compact in L 2 (m). Let Q C L 2 {m) denote the u;-limit set of this trajectory. 
As is well-known, Q is non-empty, compact, fully invariant under the evolution of (|12jl . 
and Q attracts w(t) in L 2 (m) as r — > +oo. If u) G f2, there exists a complete trajectory 
{w(r)} re R of (fT2"j) such that w(r) G Q for all r G R and w(0) = w. By Lemma l3~3l 
u)(r) = aG for all r G R, where a = J R2 d£ = J R2 wo(0 d£. Thus Q = {aG}, which 
is the desired result. □ 



3.3 Convergence in L X (R 2 ) 

We now study the behavior of the solutions of (JT2j) in L X (R 2 ) and prove Theorem II .21 and 
Proposition 11.31 

Proposition 3.4 Let w G L 1 (R 2 ) ; and let w G C°([0, oo), L 1 (R 2 )) be the solution of 
n~ty) with initial data wq. Then \w(r) — aG\i ^ as r — ► oo ; where a = J R2 u»o(£) d£. 

Proof: We know from Lemma f2. 51 that {u'(r)} r > is relatively compact in L X (R 2 ). Let 
be the cu-limit set of this trajectory. Then Q is non-empty, compact, fully invariant under 
the evolution of ([12)1 . and Q attracts w(r) in L X (R 2 ) as r — > +oo. If w G f2, there exists a 
sequence r n — > oo such that \w{r n ) — w\\ — > and «;(£, r n ) — > u>(£) for almost all £ G R 2 . 
Using (|3H|) and Lebesgue's dominated convergence theorem, we obtain 

\w{i)\ = lim M£,0| < ^l^oKe-^l 2 / 4 , £ G R 2 , (51) 

n— >oo 

since a(r n ) — > 1 as n — > oo. In particular, this shows that f2 is bounded in L 2 (m) for all 
m > 1. 

Now, since Q is invariant under the semiflow of (|12p. there exists a complete trajectory 
{w(r)} re R such that u)(r) G f2 for all r G R and w(0) = w. As we just observed, 
{w(t)} t£ r, is bounded in L 2 (m) for all m > 1. Applying Lemma 13.31 we conclude that 
u)(r) = for all r G R, where a = J„ 2 iD(^)d^ = | R2 w o(0 dC- This proves that 
Q = {aG}, which is the desired result. □ 
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Proof of Theorem OJ Let u G C°([0, oo), ^(R 2 )) be the solution of © with initial 
data lu . If we set r) = e r c<j(£ e 7 "/ 2 , e T — 1), then w G C°([0, oo), L 1 (R 2 )) is the solution 
of (|12j) with initial data u> = u . Applying Proposition 13.41 and returning to the original 
function lu(x, t), we obtain 



lim 

t— >oo 







a 



i+1 V^t+T 







which is equivalent to Q for p = 1. Next, interpolating between JSJ) for p = 1 and (JH) 
for p = oo, we obtain © for p G (1, oo). Then, using (jTUjl if 2 < g < oo or Lemma 2.1(b) 
in if q = oo, we arrive at (jlOj) . Finally, using the previous results and the integral 
equation satisfied by lu(x, t), it is straightforward to show that © also holds for p = oo. 
□ 

As we already observed, a bounded negative trajectory of (fEZJl in L : (R 2 ) need not be 
relatively compact. However, if we assume that the trajectory is relatively compact or at 
least that its a-limit set is nonempty, then we have the following result which generalizes 
Lemma 13.31 

Proposition 3.5 Let \w{r ) jvsr be a complete trajectory of Alty) in L 1 (R 2 ) ; and assume 
that w(t) has a convergent subsequence in L^R 2 ) as r — ► — oo. Then w(r) = aG for all 
r G R ; where a = J R2 w(£, 0) d£. 

Proof: By assumption, there exists w in L X (R 2 ) and a sequence r n — > — oo such that 
\w(T n ) — w\i — > and tu(£, r n ) — > for almost all £ G R 2 . Fix r G R, and take n > 
sufficiently large so that r n < r. In view of (|39jh we have for all £ G R 2 



R2 a(r - r n j 
1 

R2 a(r - r n ) 
1 



exp ( —j3 



,\£-ye 



-(r-T„)/2|2 



OT-T. 



exp I —(3 



exp ^— /5 



4a(r 



-(r-r„)/2|2 



4a(r - r n ) 

|£_ ye -(r-r„)/2|2 



u>(2/,t»)| dy 

w(y,T n ) - w(y)\dy 
- f 11 ' «•(//) <l/y 



+ 

Taking the limit n 



4a(r - r„) 

C p \w\ ie -W 2 / 4 . 

oo and using Lebesgue's dominated convergence theorem, we obtain 

\w^,T)\<C p \w\ ie -^ 2 / 4 , £ G R 2 , r G R . 

This shows that the trajectory {u>(t)} tG r, is bounded in L 2 (m) for any m > 1, and the 
result follows from Lemma f3. 31 □ 

As is clear from the change of variables ([lip, results about negative trajectories of (fT2"j) 
give information on the behavior of solutions of © as t — > 0+. In particular, Lemma l3~B1 
and Proposition 13 . 51 show that solutions of © with Dirac masses as initial data are unique 
in a certain class. A generalization of these results allows to prove Proposition 11.31 

Proof of Proposition 11.31 In view of Theorem 11.11 we can assume without loss of 
generality that T = oo. The solution u(x,t) of © can be represented as 



u>(x, t) 



F u (x,t;y,s)uj(y, s)dy , x G R 2 



t > s > . 
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Here T u is the fundamental solution of the time-dependent linear operator d t — A + u ■ 
V, and u(x,t) is the velocity field obtained from u(x,t) via the Biot-Savart law. By 
assumption, there exists K > such that < K for all t > 0. From 9j, we know 

that for any (3 G (0, 1) there exists Cp > (depending on such that 

|r„(x, t- y, s) | < i2_ exp (-/3 tlM!) , (52) 

for all x, y G R 2 and all £ > s > 0, see (J37j) . Moreover, it is shown in (see also 
Theorem 3.1 in [TH] ) that T M is a Holder continuous function of its arguments. More 
precisely, there exists 7 G (0, 1) (depending only on K) and, for any r > 0, a constant 
C > (depending only on K and r) such that that 

|r u (x, £; y, s) - r u (z, £; y', s') \<c(\y- y'V + \s - s'| 7/2 ) , (53) 

whenever t — s > t and t — s' > r. In particular, if x, y G R 2 and £ > 0, the function 
s 1 — ► r u (x, t; y, s) can be continuously extended to s = 0, and this extension (still denoted 
by r u ) satisfies (JHSJ with s = 0. 

Now, fix x G R 2 and £ > 0. Then for any s G (0, t) we have 

= / T u (x,t;y, 0)u(y, s) dy 

+ / (T u (x,t;y, s) - T u (x,t;y, 0))u(y, s) dy . 

In view of (jSHJ), the second integral in the right-hand side converges to zero as s — > 0+, 
since \u)(-, s)\i < K for all s > 0. On the other hand, since y 1— > T u (x, t; y, 0) is continuous 
and bounded, and since u(-,s) converges weakly to a5 as s — >• 0+, the first integral 
converges to aT u (x, t; 0, 0). Thus 

\u(x,t)\ = \a\\T u (x,t; 0,0) | < ^i e -^\ 2 /m , x G R 2 t>0. 

Finally, let w(£,t) = e T u(£e T/2 ,e T ) for ^ G R 2 , r G R. Then w G C7°(R, ^(R 2 )) 
is a solution of which satisfies H£,r)| < C^|a| e"^! 2 / 4 for all ^ G R 2 , r G R. 
In particular, {«;(•, r)} re R is bounded in L 2 (m) for any m > 1, hence by Lemma 13.31 
w {£> T ) = a 'G{0 for some a' G R. Clearly a' = a, and the proof is complete. □ 

Remark 3.6 As is clear from the proof, Proposition M.tA remains true if one assumes 
only that ui(-,t) stays bounded in L X (R 2 ) and that u(-,t n ) converges weakly to aS for 
some sequence t n — > 0. 

Remark 3.7 A slight extension of the techniques developped in this section allows to 
prove that the Oseen vortices are stable in the sense of Lyapunov : for any e > 0, there 
exists 5 > such that, for all initial data Wq G L x (R 2 ) with \wq — etG\\ < 5, the solution 
°f satisfies \w(r) — aG\i < e for all t > 0. Note that this does not follow from 
Proposition \3.4\ 
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3.4 Convergence rate for positive solutions 

If we restrict ourselves to nonnegative solutions of (|12jl. then combining the entropy dis- 
sipation law (J4*T|) with a few classical inequalities in information theory we can obtain an 
explicit estimate on the time needed for a solution to approach the Oseen vortex. This 
is the so-called "entropy dissipation method", which is by now a classical approach in 
kinetic theory, see for instance (2H1 d EH] ■ 

Let w e L^R 2 ), w > 0, and assume that a = J R2 w {0 d£ > 0. In information theory, 
the quantity H{w) defined in (J4l)j) is called the relative Kullback entropy of w with respect 
to the Gaussian G. Similarly, I(w) defined in ()48|) is called the relative Fisher information 
of w with respect to G. The difference between the entropy H(w) and its minimal value 
H(aG) is bounded from below by the Csiszar-Kullback inequality 

— \\w-aG\\ 2 Ll < H(w)-H(aG) , (54) 

and from above by the Stam-Gross logarithmic Sobolev inequality 

H(w) - H(aG) < I(w) . (55) 

Assume now that w € L 2 {m) PI S + for some m > 3, and that a = J R2 ^o(0 d£ > 0. 
Let w G C°([0, oo), L 2 {m)) be the solution of (fT^j) with initial data w . Combining (|4Tjl 
with we immediately obtain 

H(w(t)) - H(aG) < (H(w ) - H(aG)) e' T , r > . 

Applying (jo^j) . we conclude that 

\\w(t) -aG\\ L i < V2a~{H(w ) - H{aG)) 1/2 e~ r/2 , r>0. (56) 

This shows that w{r) converges to aG at the rate e~ T / 2 , which is optimal in general (see 
Section EJ). Moreover, (|5T)|l gives an explicit estimate of the prefactor in terms of the 
initial data. In particular, this provides an explicit upper bound of the time needed for 
the solution to enter a given neighborhood of the vortex. 

Unfortunately, we do not know how to extend the entropy dissipation method to the 
general case where the vorticity may change sign. In the next section, we obtain local 
convergence rates by studying (J 12)) in a neighborhood of the family of Oseen vortices, but 
this approach does not provide any explicit estimate in the sense of (|56|). 

4 Local convergence rates 

From the results of the previous section we know that any solution of the Navier-Stokes 
equation whose initial vorticity distribution lies in L X (R 2 ) will converge toward the Oseen 
vortex with the same total vorticity. In the present section we show that for solutions 
in the weighted space L 2 (m) with m > 1 we can derive precise estimates on the rate at 
which solutions approach the vortices. 

Our analysis proceeds by first analyzing the linearization of fT2|) at a vortex solution. 
We prove estimates on the location of the spectrum of linearized operator which in partic- 
ular imply that the vortex solutions are spectrally stable for all values of the circulation 
Reynolds number. We then show that these bounds also imply decay estimates for the 
full nonlinear evolution in a neighborhood of the vortex. 
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4.1 Eigenvalue estimates 

Fixing a G R and linearizing (|T2j) around w = qG, we find 

d t w + ct(v G • Vw + v ■ VG) = Aw + -£ ■ Vw + w , (57) 

where as usual v is the velocity field associated to w via (fUj) and v G is the velocity field 
of the Oseen vortex. This equation can be rewritten as d t w = £w — a Aw, where 

Cw = Aw + -£ ■ Vw + w , Aw = v G • Vu> + v ■ VG , (58) 

see (J2H1), (El- 

The linear operator £ in L 2 (m) is studied in detail in Q14J, Appendix A). It is defined 
on the maximal domain 

V m {L) = |w G L 2 (m) Amj + ^- Vmj6 L 2 (m)| . 

If w G T> m (£), one can show that Aw G L 2 (m), so that T> m (C) C H 2 (m). The essential 
spectrum of £ is given by 

1—171' 



<C{£) = {\ec\Be{\)<—^} 



In addition, C has a sequence of eigenvalues 0, —1/2, —1, . . . whose eigenfunctions are 
rapidly decreasing at infinity. 

Because of the spatial decay of v G and G, the operator A is a relatively compact 
perturbation of £ and hence d£f(£ — aA) does not depend on a. In particular, we can 
always push this essential spectrum far away from the imaginary axis by taking m > 
sufficiently large. Thus the spectral stability of the vortex solutions will be determined 
by the isolated eigenvalues of C — a A in L 2 (m). As we shall see, the corresponding 
eigenfunctions have a Gaussian decay at infinity so that, in contrast to the essential 
spectrum, these isolated eigenvalues do not depend on m. 

We next observe that, due to symmetries of equation (J2J), some eigenvalues of C — aA 
are in fact independent of a. For instance, if m > 1, then A = is a simple eigenvalue of 
£ in L 2 (m), with eigenfunction G. Since v G • VG = 0, it is clear from (|58|) that AG = 0, 
so that is an eigenvalue of £ — a A for any a G R. This zero eigenvalue is due to the fact 
that the Oseen vortices form a one-parameter family of equilibria of (|12|). The associated 
spectral projection P reads (for any a) 



JR 2 



Thus, it will be sufficient to study the spectrum of £ — aA in the spectral subspace L^m) 
defined in (|15j) . which by (0) is also invariant under the nonlinear evolution (|28j). 

Similarly, if m > 2, A = —1/2 is a double eigenvalue of £ with eigenfunctions Fi,F 2 , 
where Fj = —djG (J = 1,2). Differentiating the identity v G • VG = with respect to £j, 
we see that AFj = (j = 1, 2). It follows that —1/2 is still an eigenvalue of £ — aA for 



22 



any a G R. This eigenvalue originates in the translation invariance of (J2J) with respect to 
£ G R 2 . The associated spectral projection Pi reads (for any a) 

(p lW )(o = f 1 (o [ exo d^ + f 2 (o / exo . 

JR 2 JR 2 

Thus, it is again sufficient to study the spectrum of C — aA in the spectral subspace L 2 (m) 
defined by 

L\(m) = [we L 2 (m) | jf^ £^(0 d£ = for j = 1, 2} , (59) 

which by (jUJ) is also invariant under the nonlinear evolution (|28|) . 

Finally, if m > 3, A = — 1 is a triple eigenvalue of £ with eigenfunctions AG, (df—d 2 )G, 
and d\diG. Since AG = ^(|£| 2 — 4)G is radially symmetric, it is clear that A(AG) = 0, 
so that —1 is still an eigenvalue of C — a A for any a G R. This is due the fact that (J2J) 
is autonomous and invariant under the rescaling u(x,t) \— > X 2 u(Xx, X 2 t). However, as we 
shall see, the eigenvalue —1 is simple if a 7^ 0. The associated spectral projection P 2 reads 
(for any q^O) 

(P 2W )(0=AG(£) / V| 2 - 4MO d£'. 

Thus, if a 7^ 0, it is sufficient to study the spectrum of C — aA in the spectral subspace 
L|(m) defined by 



Z|(m) = { W G L 2 (m) | J \tfw{i) d£ = o} , 



(60) 



which (as can be verified by a direct calculation) is also invariant under the nonlinear 
evolution 

The principal result of this subsection is: 

Proposition 4.1 Fzx m > 1 and a G R. T/ien any eigenvalue X of C — aA in Lo(m) 
satisfies 

_^ / ^ \ / 1 1 — m\ , . 

Re (A) < max — , . (61) 



2 ' 2 

If moreover m > 2, then any eigenvalue X of £ — aA in L 2 (m) satisfies 

Re (A) < max(-l , — — J . (62) 

Remark 4.2 In view of the preceding remarks estimates \61)) and \6fy) are sharp. If 
m > 3 and a 7^ 0, the proof shows that any eigenvalue of C — a A in L\{m) satisfies 
Re (A) < —1, but we are not able to give a sharp estimate in that case. Numerical 
calculations in \25$ indicate that the eigenvalues that are not frozen by symmetries have a 
real part that converges to —00 as \a\ —* 00, thereby suggesting that a fast rotation has a 
stabilizing effect on the vortex. Proposition \4-l\ shows rigorously that at the spectral level 
perturbations of the vortex solutions decay at least as fast when a is large as when a = 0. 
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To prove Proposition 14. 11 we proceed in three steps. First, we observe that the linear 
operators C and A are invariant under the group of rotations SO (2). Thus, using polar 
coordinates in R 2 and expanding the angular variable in Fourier series, we reduce the 
eigenvalue equation for the operator C — aA to a (nonlocal) ordinary differential equation 
in the radial variable. Next, a careful study of this differential equation reveals that, 
if A G C is an isolated eigenvalue of £ — aA, the corresponding eigenfunction has a 
Gaussian decay at infinity. Finally, we prove that the operators C and A are respectively 
self-adjoint and skew-symmetric in a weighted L 2 space with appropriate Gaussian weight, 
and Proposition 14. II then follows from elementary considerations. 



4.1.1 Polar coordinates 

Fix m > 0. For any n G Z, let P n be the orthogonal projection in L 2 (m) defined by 

(P n w)(r cos 8, r sin 8) = uo n (r)e vnd , where 

u) n (r) = — [ w{r cos 8, r sin 8)e~ ine dd . (63) 
2vr J 

Clearly, P n P n > = S nn 'P n and J2nez Pn = 1- ^ w e L 2 (m), the functions u> n : R + — ► C in 
(jnHl) belong to the Hilbert space 



Z(m) = \uj : R, -> C 



j r(l +r 2 ) m \uj(r)\ 2 dr < ooj 



(64) 



For any n G Z, let £ n be the linear operator on Z(m) defined by 

C n uj = d 2 u+ (^ + ~y r u+ ■ (65) 

Let also A n be the (bounded) linear operator on Z(m) defined by Ao = and 

A n u = in(ipuj — gQ) , n^O, (66) 

where 
and 

1 / fWz\H , , /' 0C /r\N 



n < r » = U W zwWdz+ / W 2 " (2)d ^ • (67) 

It is easy to see that the operator A n is indeed well-defined: 

Lemma 4.3 If n G Z* and u G Z(m) for some m > 0, then J6'7| ) defines a continuous 
function Q : R + — > C. Moreover, Q(r)/r converges to zero as r — > oo and zs oi mosi 
logarithmically divergent as r — > 0. 

Proof: The proof is straightforward using (j64jl . (|6Tjl and Holder's inequality. □ 
We now show that £ n and A n are the expressions of £ and A in polar coordinates: 
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Lemma 4.4 The operators C and A commute with the projections P n . If n G Z and 
w G D m (C) for some m > 0, then 

(CP n w)(r cos 0,rmn0) = e ine (£ n cj n )(r) , (68) 
(AP„w)(rcos#,rsinfl) = e ine (A„cj n )(r) , (69) 

where u n (r) = e~ m0 (P n w)(r cos9,r sin 9). 

Proof: All we need is to prove (|5Sjl and ()69|) . The first relation follows from (J55J) by an 
elementary calculation. To prove (J69j) . assume that w(r cos 9, r sin 6) = u n (r)e for some 
nGZ and some u n G Z(m). Then 

v G ■ = - e- r2/4 )^ u; = ine in6 cpu n . 

lur r 

On the other hand, the velocity field v corresponding to w satisfies div 2 — d 2 v\ = w, 
diVi + d 2 v 2 = 0. In polar coordinates, these relations become 

-dJrvg) - -dgv r = w , -d r (rv r ) + -dgvg = . 

We look for a solution of the form v r = v r (r) e me , vg = vg(r) e me . Then 

(rvg)' — inv r = roj n , (rv r )' + invg = . 
Eliminating vq, we find the following ODE for h = rv r : 

(rh'Y h + inruj n = . 

r 

The general solution is: 

k{r) = 2H \J Vr) ZUJn ^ dz + j {-) ^n(z)dzj+A 1 r n + A 2 r- n , 

where A\,A 2 G C. Since we want a velocity v r = h/r that is locally integrable and 
converges to zero at infinity, we must choose Ai = A 2 = 0. Setting Q = -^h = h^rv r , we 
finally obtain: 

v ■ VG = -~rv r g = -ine ine gVt . 
This concludes the proof of ()69j). □ 

4.1.2 Gaussian decay of eigenfunctions 

The aim of this paragraph is to prove: 

Lemma 4.5 Fix m > 0, and assume that w G T> m (C) satisfies (C — aA)w = fiw, where 
a G R and Re (fi) > ^j— . Then there exist C > and 7 > st/c/j i/iai 



K0l<C(i + |£| 2 ) 7 e-i«' 2 / 4 , eeR 
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Proof: We use the notations of the preceding paragraph. According to Lemma f4. 41 we 
can assume that the eigenfunction w satisfies w = P n w for some nGZ. Thus, there exists 
u) G Z(m) such that w(r cos 8, r sin 6) = u(r)e md . In view of (|65|). and Lemma f4.4t u 
satisfies the (inhomogeneous) ordinary differential equation 



/ r 






i) J(r) + ( 







rr 



(70) 



where Q is defined in (|57j) . 

The basic idea is now to use the classical results of Coddington and Levinson [TTJ to 
show that, for r large, any solution of (J7U|) either decays like u(r) ~ r 27 e _r / 4 for some 
7 > 0, or like u{r) ~ r 2At_2 in which case it cannot belong to L 2 (m). However, a certain 
amount of preliminary work is needed in order to bring (j70j) into a form to which we can 
apply the results of JTTj. 

We begin by introducing new variables /, F and t via the definitions 

oj( r ) = f(r 2 /A) , n(r) = F(r 2 /4) , t = r 2 /4 . 
In terms of these new variables, (|70p takes the form 

f"(t) + (l + /'(t) + - a(t)\ f(t) + b(t) = , t > , (71) 



where 



2 



Recall that we are interested in the behavior of solutions of this equation for t large. 
We first consider the behavior of the homogeneous part of this equation and then construct 
the solution of the full equation via the method of variation of parameters. 

Lemma 4.6 The linear, homogeneous equation 

f(t) + (i + 7) /'(*) + {^r - a (t)) m = - 

has two linearly independent solutions <fi(t) and (pzfy) such that 



t— >oo 



<Pl'(t)J \0J ' t^oo \<p 2 '(t)J \-l 



Proof: If we define x(t) = (^p^j we can rewrite the differential equation in the lemma 

as 

x'(t) = (A + V(t) + R(t))x(t) , 

where 



t t 



a(t) 



But in this form, the lemma follows immediately from (jTTj, Theorem 3.8.1, p. 92). □ 
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We now derive the asymptotic form of the solution f(t) of (|71|) by applying the method 
of variation of parameters. Set 



m = A(t)<p 1 (t) + B(t)<pi(t) , (72) 
where A'(t)ipx(t) + B'(t)(p 2 (t) = 0. Then differentiating / and using (fTTj) we find 

A\t)<f/ 1 (t) + B'(t)<p' 2 {t) + b{t) = 0. 
Solving for A' and B' we obtain 

'A'(t)\ _ 1 / b{t)ip 2 {t) 



B'(t)J W(t) V-6(t)^(t) 

where W(t) = — fe~* is the Wronskian determinant of (fx and tp 2 . Integrating both sides 
of this equation we find 



A(t) 



A x - J se s b(s)ip 2 (s) ds , B(t) = B x + J se s b(s)ipx(s) ds . 



Recall from the definition of b(t) that e s b(s) = *is£l£2. From Lemma 14.31 we know that 
F(s)/s = n(2s/s)/s converges to zero as s — > oo, so that e s b(s) is bounded for s > 1. 
Thus, the asymptotic behavior of (p 2 implies that Ax(t) — > A^ as t — > oo. If 7^ 0, then 
f(t) ~ Aoot' 1-1 as t —>■ 00 and hence (reverting to the original polar coordinates) u{r) ~ 
Woo"?" 2 ^ 2 - But since Re (//) > ^p, this would imply that / °°r(l + r 2 ) m \u(r)\ 2 dr = 00 
and this in turn would violate the hypothesis that the eigenfunction w is in L 2 (m). Thus, 
Aoq = 0, and A(t) = f°° se s b(s)ip 2 (s) ds from which we conclude that \A(t)\ < Ce'H' 1 , for 
some 7 > 0. In analogous fashion one proves that |-B(i)| < CP . Inserting these bounds 
on A and B into ()72j) and using the asymptotic estimates on ipx and ip 2 we conclude that 
there exists 7 > such that \f(t)\ < Ct 7 e~*, for t > 1, or 



UJ 



(r)| < Cr 27 e- r2/4 , r>l 



This is the desired estimate, since |w(£)| = 1^(1^1)1- The proof of Lemma 14.51 is now 
complete. □ 

4.1.3 Localization of eigenvalues 

Let X denote the (complex) Hilbert space 

X={we L 2 (R 2 ) I G~ 1/2 w e L 2 (R 2 )} , 
equipped with the scalar product 

(wi,w 2 ) x = I -^7-rWi(£)w 2 (Z)d£ . 

JYi? <~H<U 

We also introduce the closed subspaces X ,X 1 defined by 

x = {wex\ y ^(e) de = 0} = XHL 2 (m) , 

Xx = {wGl J Zjwffl d£ = for j = 1, 2} = X n L\{m) . 
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According to Lemma l4.5[ if /i is an eigenvalue of £ — a A in L 2 (m) with Re (/x) > , . 
the corresponding eigenfunction belongs to X. This result is very useful because both 
operators £ and A have nice properties in this space. 

Lemma 4.7 The linear operator £ is self-adjoint in X, and £ < 0. Morevoer, £ < —1/2 
on X and £ < — 1 on X\ . 

Proof: Define L : V(L) L 2 (R 2 ) by V{L) = {ip E H 2 (R 2 ) \ \£\ 2 i/j E L 2 (R 2 )} and 

L = G- l ' 2 £G 1 ' 2 = A- S + 1 . 

16 2 

In quantum mechanics, the operator — L is (up to numerical constants) the Hamiltonian 
of the harmonic oscillator in R 2 . As is well-known (see for example [E])> L is self-adjoint 
in L 2 (R 2 ) and o~(L) = {—n/2 \ n E N}. By construction, the operator £ : V{£) — > X 
with domain 

V{£) = E X \£\ 2 w E X , Aw + -£ • Vw E x| , 

is thus self-adjoint in X with the same spectrum. In particular, £ < 0. Now, observe that 
is a simple eigenvalue of £ with eigenfunction G, and that X is just the orthogonal 
complement of the eigenspace RG in X. Thus Xq is stable under £ and the restriction 
of £ to Xo is a self-adjoint operator satisfying £ < —1/2. Similarly, one can show that 
£ < -1 on Xi. □ 

Lemma 4.8 The linear operator A is skew- symmetric in X. 

Proof: Since AG = and since the subspace Xo is stable under A, it is sufficient to 
show that A is skew-symmetric on X . Let w,w E X fl T>(£), and denote by v, v the 
corresponding velocity fields. Without loss of generality, we assume that w, w (hence also 
v, v) are real functions. Then 

(w, Aw)x = f 2 (J=< ™v G ■ Vw - ^w(v • £)J d£ , 

because VG = — Observe that G _1 v G is a divergence free vector field, so that 

1 f 1 

— w\ G ■ Vw d£ = - / — w\ G -Vwd£ . (73) 

R 2 G J R 2 G 

On the other hand, the following identity is easy to check: 

w(v ■ + w(v ■ = (^dt - 6^2)K^2 + v 2 v 1 ) + (^d 2 + &di){v 2 v 2 - fi^i) . 

Since w,w have zero mean, it follows from (|14j. Appendix B) that v,v decay at least 
like 1/|£| 2 as C, — > oo. Thus, integrating both sides, we obtain 

(w(v-0+w(v-0)d^ = . (74) 

R2 

Combining (J73J) and (f74"|) . we see that (w, Aw)x + (Aw,w)x = 0. □ 
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Proof of Proposition 14.11 Fix m > 1 and assume that A is an eigenvalue of £ — aA in 
Lg(m) with Re (A) > By Lemma l4.5( there exists a nonzero w G X fl T>(£.) such 

that (£ — «A)u> = Aw. In particular, 

X(w,w) x = (w,Cw) x - a(w,Aw) x , 

hence 

Re(A)(iy,w)x = (w,£w) x < --(w,w) x , 

since A is skew-symmetric and £ < — 1/2 on X . Thus, Re (A) < —1/2. 

Similarly, if m > 2 and A is an eigenvalue of £ — aA in L\{m) with Re (A) > 
there exists a nonzero w G X] fl £*(£) such that (£ — aA)w = Aw. Proceeding as above 
and using the fact that £ < — 1 on Xi, we obtain Re (A) < — 1. □ 

Remark 4.9 If m > 3 and A is an eigenvalue of £ — aA m -L^ 771 ) Re (A) > , 
there exists a nonzero w G P(£)nL|(m) suc/t t/iat (£ — aA)w = Aw. The argument above 
shows that Re (A) < — 1, and that Re (A) = —1 if and only if Cw = —w. But this implies 
that w is a linear combination of (d\ — and did 2 G, and a direct calculation shows 

that no such w can be an eigenfunction of £ — a A i/a^O. Thus any eigenfunction A of 
£ — a A in L\(rn) satisfies Re (A) < —1 if a ^ 0. 



4.2 Bounds on the linear evolution 

In this subsection we prove that the eigenvalue estimates of the previous subsection imply 
analogous bounds on the linear evolution. Fix a G R and consider the linear equation 
d T R = CR — aAR which is the linearization of (fT^j) a the vortex aG, see fl2~%j) . (JSTj) . The 
corresponding integral equation reads 

R(r) = S(t)R -a [ S{t- s)AR(s) ds , (75) 
Jo 

where <S(t) = exp(r£). Proceeding as in Section |2~H it is straightforward to show by a 
contraction mapping argument that this equation defines a strongly continuous semigroup 
T a (r) in L 2 (m) for any m > 1, namely R{t) = 7^(r)i? - We first prove that T a (r) is a 
compact perturbation of S(t). 

Lemma 4.10 Letm > 1. The linear operator JC a (r) = T a (r) —S(t) is compact in L 2 {m) 
for each r > 0. 

Proof: All the necessary estimates are already contained in Section |2~T1 Observe that the 
term K, a (r)R in (|73j) is precisely what we called R 2 (r) + R 3 (t) in the proof of Lemma l2~2l 
Repeating the estimates proved there, we obtain 

\\}C a (T)R \\ m+l < C(T)\\Ro\\ m , r>0, 

for some C(r) > 0. Similarly, the result of Lemma 12.11 applies to T a (r)Ro, hence to 
JC a (r)Ro, and yields 

\\VIC a {T)R \\ m <-^-\\R \\ m , r>0, 
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where a(r) = 1 — e r . The conclusion then follows from Rellich's criterion. □ 

We now use this result to bound the essential spectrum of T a (r). We recall that A G C 
is in the essential spectrum of a linear operator L if A is not a normal point for L, i.e. if A 
is not in the resolvent set of L and is not an isolated eigenvalue of L of finite multiplicity, 
see e.g. [20] • For any r > 0, we denote by £>(r) the closed disk of radius r centered at the 
origin in C: 

B(r) = {z G C | |*| < r} . 

Lemma 4.11 Let m > 1. For any r > 0, the essential spectrum of T a {r) in L 2 (m) 
satisfies 

<CP*«(r)) = B(e^- m » 2 ) . (76) 

Proof: Fix r > 0. The results of Appendix A) imply that the essential spectrum 
of <S(t) in L 2 (m) is exactly ^(e^ 1 -" 1 )/ 2 ). Since T a (r) is a compact perturbation of S(r), 
it follows from ([20J, Theorem A.l) that the complement of the B(e T ^ 1 ~ rn ^' 2 ) in C either 
consists entirely of eigenvalues of T a (r), or entirely of normal points for T Q (r). In our 
case, the first possibility is excluded. Indeed, assume that A G C is an eigenvalue of 
T a (r) with |A| > 1. By the spectral mapping theorem, A = e v where v is an eigenvalue of 
C — aA in L 2 {m) with Re {y) > 0, which contradicts Proposition 14.11 Thus, cr^ s (7^(r)) C 
j6(e r ( 1 ~ m )/ 2 ), and since S(r) is also a compact perturbation of T a (r) the same argument 
shows that a^(T a {r)) = B{e T ^- m ^ 2 ). □ 

By construction, the spectral subspaces L^{m) and L\(m) are left invariant by the 
semigroup T a (r). Combining Lemma 14.111 with the eigenvalue estimates of Section f4.1l 
we obtain precise bounds on the growth of T a (r) in these subspaces: 

Proposition 4.12 Fix a G R. Assume that either 

a) 0</i<l/2 and Rq G Lg(m) for some m > 1 + 2fi, or 

b) 1/2 < fi < 1 and Rq G L\(m) for some m > 1 + 2/x. 
T/ien there exists C > (independent of Rq) such that 

\\r a {r)Ro\\ m < Ce-^\\R \\ m , r > . (77) 

Proof: Assume first that < fi < 1/2. By Lemma f4. Ill the essential spectrum of T a (l) 
in L 2 {m) satisfies < s (T a (l)) = B(e^' m ^ 2 ), and e^ 1 "^/ 2 < e"" since m > 1 + 2/i. If 
A G <x(T a (l)) satisfies |A| > e^ 1 -" 1 )/ 2 , then A is an isolated eigenvalue of T a (l) and (by 
the spectral mapping theorem) there exists an eigenvalue v of C — a A in Lp(m) such that 
e u = A. Applying Proposition 14.11 we obtain Re (y) < 1/2, hence |A| < e^ 1 / 2 < e~ M . Thus 
the spectral radius of T Q (1) in Lo(m) is strictly less than e _At , and (|77j) follows (see e.g. 
[T3] . Proposition IV. 2. 2). A similar argument proves (J77|) if Rq G L 2 (m) and 1/2 < // < 1. 

Now, assume that fj, — 1/2 and m > 2. Any i? G L^m) can be decomposed as 
#0 = /^Fx + /3 2 F 2 + .Ro, where ^ = J R2 df , Fj = -djG, and R G L 2 (m). It follows 
that 

T a (r)R = e^l 2 (p x F x + (3 2 F 2 ) + T a (r)R , 

and we already know that ||^,,(r)i?o ||m < Ce~ yr ||i?o||m for some v > 1/2. Thus (J77J) 
holds and is sharp in this case. A similar argument shows that f)77j) holds with /x = 1 if 
i?o £ L\{m) for some m > 3. □ 



30 



Finally, we will need in the following subsection L p -L q estimates of T a { T ) and its 
derivatives, in the spirit of ()23j). 



Proposition 4.13 Under the assumptions of Proposition "J^T^ if Rq satisfies in addition 
b m R G L q (R 2 ) for some q G (1,2), then 

\\r a (r)Ro\\ m < ° e , \\VT a (T)R \\ m < \b m R \ q , r > . 

a(r)i 2 a(r)i 

Proof: Fix 1 < q < 2 and assume that b m Ro G L q (R 2 ). Given T > 0, we consider the 
function space 

Y = {R G C°((0,T],H 1 (m)) \ \\R\\ Y < oo} , where 
\\R\\ Y = snp \a(T)^\\R(r)\\ m + a(r^\\VR(r)\\ m } . 

0<t<T *> J 

Using (j2SJ) and proceeding as in the proof of Lemma l2~T| it is straightforward to show that, 
if T > is sufficiently small, there exists C > such that (|73j) has a unique solution in the 
ball of radius C\b m Ro\ q centered at the origin in Y. This proves the desired estimates for 
< t < T, and the general case follows from Proposition 14. 12l if one uses the semigroup 
property and the fact that the subspaces L^m) and L^(m) are left invariant by T a {r). □ 

4.3 Bounds on the nonlinear evolution 

In this subsection, we show that our estimates on the linear semigroup T q (t) generated 
by C — a A imply similar bounds on the full nonlinear evolution (|28j). As is easy to verify, 
the integral equation (|3U[) satisfied by R(t) is equivalent to 

R(r)=T a {r)Ro- I T a (r - s)(v R (s) ■ VR(s)) ds . (78) 
Jo 

The following result implies Proposition 11.61 as a particular case: 



Proposition 4.14 Under the assumptions of Proposition 4-l 6 A the solution R(t) of 
with initial data Rq satisfies \\R(r)\\ m = 0(e~^ T ) as r — > oo. 

Proof: Let R G C°([0, +oo), Ll(m)) be the solution of (J2H) with initial data Rq. We know 
from Proposition 11.51 that ||i?(r)|| m converges to zero as r — » oo, hence we can assume 
without loss of generality that ||i2o|| m is small. Given T > 0, we define 

M(T) = max{ sup e^\\R(T)\\ m , sup a(r)2 e^ r || VR(r)\\ m } . 

0<r<T 0<t<T J 

By Lemma O M(T) < oo and there exists C > such that M(T) < C \\R \\ m if T 
is sufficiently small. Applying Propositions 14.121 and 14.131 and using the fact that the 
subspaces Li^{m) and L\(m) are stable under the nonlinearity • V-R, we obtain from 

||i*(r)|| m <Ce-' M -||i2o|| m +C f e'^^ l —--\b m v R { s ) ■ VR(s)\ q ds , 

Jo a(r—s)i 2 
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where 1 < q < 2. As in the proof of Lemma \2. 11 we can bound 

\b m v R (s) ■ VR(s)\ q < C\\R(s)\\ m \\VR(s)\\ m < —-^M(Tf 

a(s)2 



hence 



\\R(r)\\ m < e-^(c\\Ro\\ m +CM(T) 2 / — r ds 

v Jo a(r—s)i 2 a(s)a 

< e^CipolU + C 2 M(T) 2 ) , 0<r<T, 



where C\,C 2 > are independent of T. Without loss of generality, we assume in what 
follows that C\ > Co- Differentiating (fTSj) and using similar estimates, we also obtain 



p H" / \ 

\\VR(T)\\ m <—- T (C 1 \\R Q \\ m + C 2 M(T) 2 ) , 0<r<T. 
a(rj2 v ' 

Summarizing, we have shown 

M(T) < CiH^IU +C 2 M(T) 2 . (79) 

Now, assume that Rq is small enough so that 4CiC2||i?o|| m < 1- As C\ > Cq, we also 
have Co||-Ro||m < M, where 



M -- 

2C 2 

Thus M(T) < M for T > sufficiently small, and since M(T) depends continuously on 
T it follows from (O that M(T) < M for all T > 0. In particular, ||i?(r)|| m < Me""" 
for all r > 0. □ 

Remark 4.15 Similarly, if Rq G L^m) for some m > 3, the solution R(t) of \2§) with 
initial data Rq satisfies \\R(r)\\ m = 0(e~^ T ) as r — ► oo, for some fi > 1 depending on a. 
We know from Remark\4.9{ that /i> 1 «/«/0, but we have no sharp result in that case. 



Remark 4.16 It follows directly from the proof that the Oseen vortices are stable equi- 
libria in L 2 {m), see Remark \3~l\ 

We conclude this section by showing that, in the case where the total vorticity is 
nonzero, there is no loss of generality in assuming that the perturbations of the vortex have 
vanishing first order moments. Indeed, assume that 1/2 < // < 1 and that Wq G L 2 (m) 
for some m > 1 + 2/i. Let w G C°([0, oo), L 2 (m) be the solution of (jl2j) with initial 
data Wq. For any b G R 2 , the function w(£,t) = w(£ + &e -r / 2 ,r) is again a solution 
of (|12|) (because the original equation (J2J) is translation invariant in x G R 2 ). If a = 
| R 2 %(0 ^ 7^ 0, we can choose {b u b 2 ) = a~ 1 {f3 1 , j3 2 ), where $ = J R2 6^o(0 Then 
J R 2 £,iw(£,, 0) d£ = for i — 1,2, so that u>(-, 0) — aG G L\{m). Applying Proposition 14. 141 
we obtain ||w(t) — a(j|| m = 0(e~^ T ) as r — ► oo. Returning to the original function ■?/;(£, r) 
and using a straightforward Taylor expansion, we obtain the second order asymptotics 

\\ w (r) - aG - (ftFi + /3 2 F 2 ) e- r / 2 || m = 0(e~^), r -> oo , (80) 

where i*} = —djG, j = 1, 2. As was already mentioned, this result still holds when a = 
except that, if \i = 1 and ((3i, fl 2 ) ^ (0,0), the right-hand side of (JSUj) should be replaced 
by 0(re- T ). 
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Remark 4.17 Assume that m > 3, a ^ 0, and J R2 |£| 2 Wo(£)d£ ^ 0- V we translate and 
rescale w appropriately, we can produce a new initial condition w with w — aG G L\(m). 
By Remark \4-15\ the corresponding solution then satisfies \\w(t) — aG\\ m = 0(e~ >1T ) as 
t — > oo for some \i > 1. Moreover u>(£, r) and w(£,t) are linked by a simple relation, 
due to the fact that the original equation HJ) is translation and dilation invariant. Using 
this relation, we obtain that the the next correction to the asymptotics \8U\) is of the form 
j (AG) e~ T , for some 7 G R. 

Remark 4.18 The connection of the translation invariance of the Navier-Stokes equation 
with the decay associated to the first moment of the vorticity seems first to have been 
remarked upon by Bernoff and Lingevitch, J3J^. The connection between symmetries of 
the linear and nonlinear heat equation and Burgers' equation and the decay rates of the 
long-time asymptotics of solutions of these equations was systematically explored in VJl^ 
and JZEfl. 
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